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Abstract 

We study the quasiparticle transport coefficients in disordered d-w&ve super- 
conductors. We find that spin and charge excitations are generally localized 
unless magnetic impurities are present. If the system is close to a nesting point 
in the impurity-scattering unitary limit, the tendency towards localization is 
reduced while the quasiparticle density of states gets enhanced by disorder. 
We also show that the residual repulsive interaction among quasiparticles has 

a delocalizing effect and increases the density of states. 
72.15.Rn,74.20.-z,74.72.-h,74.25.Fy 
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The discovery of ci-wave superconductivity in doped cuprate materials has stimulated a 
substantial experimental and theoretical activity focused on the role played by disorder in 
superconductors with gapless quasiparticle excitations. 

Indeed this type of superconductors raises several questions which are worth being ad- 
dressed. For instance, charge is not anymore a conserved quantity in the superconducting 
phase, contrary to spin and energy. Therefore, while quasiparticle spin and thermal den- 
sity long wavelength fluctuations are diffusive in the presence of disorder, charge density 
fluctuations do not diffuse, realizing a sort of spin-charge separation in the hydrodynamic 
limit. As it has been shown in Ref. [|TJ, the effective field theory, which describes the spin 
and thermal diffusive modes, has still the form of a non linear a-model, like in the standard 
Anderson localization. The /3-function of the appropriate cx-model implies that in two di- 
mensions the spin diffusion constant vanishes. However, in view of the previously discussed 
independence between charge and spin long-wavelength excitations, spin and charge might 
behave differently with respect to localization, resulting into a quite intriguing scenario. 

In this work we rexamine the role of quantum interference corrections to the transport 
properties of quasiparticles in a d-wave superconductor, including besides spin also quasipar- 
ticle charge conductivity. Moreover, to keep our analysis as general as possible, we consider 
several possible symmetries of the underlying model. Namely, we study models with and 
without time reversal invariance, as well as in the presence or absence of the so-called sublat- 
tice symmetry |fj], which might be relevant close to a nesting point in the unitary scattering 
limit. In all the above symmetry classes, we derive the renormalization group (RG) equa- 
tions for the charge and spin conductivities as well as for the density of states (DOS). We 
find that the spin and charge conductivities behave similarly in the absence of sublattice 
symmetry, and differently otherwise, in which case the spin conductivity stays finite while 
the charge one vanishes. The impurity induced DOS in the absence of sublattice symmetry 
is suppressed for potential scattering, in agreement with Ref. ||, but enhanced in the pres- 
ence of spin-flip scattering, which also gives a delocalizing correction to conductivities. On 
the contrary, if sublattice symmetry holds, the DOS diverges. Finally, we also consider the 

2 



quasiparticle interaction. Consistently with the charge not being a conserved quantity, we 
find that the singlet particle-hole channel does not contribute, while the triplet particle-hole 
and the Cooper channels have a delocalizing effect, which also enhances the DOS. 
We consider the following Hamiltonian on a square lattice 

{ij) a 

+E[ A ij(4Ai+ c }A)+H.c], (i) 

{ij) 

where (ij) means that the sum is restricted to nearest neighbor sites, c\ a creates an electron 
with spin a =f , \ at site i, while c ia annihilates it. We take a gap function Ajj of ci-wave 
symmetry. The hopping matrix elements are independent random variables and satisfy 
tij = tji G Tie, and 4>ij — ~4 > jii with 0^ zero or finite depending whether or not time 
reversal invariance holds. The spectrum of the Hamiltonian (|I]) posseses a nesting property. 
Namely, given a generic eigenfunction with energy E and amplitude 4>{i)e at site % = (n, m), 
the operator 0^{%)p = ( — l) ri+rra <^>(i) which shifts by (tt, 7r) the momentum, generates the 
eigenfunction with energy — E. This implies an additional symmetry (chiral symmetry) at 
E = 0. The localization properties are quite different whether chiral symmetry holds, which 
corresponds to the Fermi energy Ep = (half- filling) , or broken, Ep ^ 0. In the latter case, 
the localization properties of (|J) are analogous to models in which on-site disorder is present 
or next-nearest neighbor hopping is included, which break chiral symmetry everywhere in 
the spectrum. 

We analyse the disordered Hamiltonian (|I|) by using the replica trick method within the 
path integral formalism ||. We introduce the vector Grassmann variables q and q with 
components (k,a,p,a an d Ci,a,p,ai where i refers to a lattice site, a to the spin, p = ± is the 
index of positive (+o>) and negative (—u) frequency components, and a = 1, . . . , n is the 
replica index, as well as the Nambu spinors 



and = [C^J/j]*, with the charge conjugation matrix C = io y T\. Here and in the following, 
the Pauli matrices cr& (b = x,y,z) act on the spin components, s& (b = 1,2,3) on the 
frequency (retarded/advanced) components, and T& (b = 1, 2,3) on the Nambu components 
c and c. The action corresponding to (|l|) is 

S = E *i ( %e : + iAij-raS! - i^-wsg) ^ (2) 

ij 

As in the standard Abrikosov and Gorkov approach to superconductivity, the gap function 
couples fermions with opposite frequency u>. If magnetic impurities are present, we must 
add to @ an additional spin-flip scattering term J2i ^iTsa ■ S^i, being Si the impurity 
spin. The same term, with Si = B, gives the Zeeman splitting in the presence of a constant 
magnetic field, which also breaks time-reversal invariance. From the action (0) we can 
deduce the proper symmetries of the field theory. We consider two different global unitary 
symmetry transformations, one for sublattice A and another for B, namely ^a ~^ Ta^a 
and — > Tb$ b . For frequency u = 0, the action is invariant if CT a C 1 HabTb = Bab- 
By writing T4 = exp w °+ w,i ; Tb = exp Wo ~ W3 ; with antihermitean Ws, we find that 



r 2 si,W (3) = and r 3 ,iy (3) 



0, the 



CW l C l = -W and CW^C* = W 3 , as well as 
latter condition valid if time reversal invariance is broken. If magnetic impurities are present, 
then [r 3 cx, W ] = {r 3 a, W 3 } = 0, while r 3 a • B, W = |r 3 <? • B, W 3 | = in the presence of 
a constant magnetic field. Since the frequency u acts as a symmetry breaking field, we have 
to further impose that {Wo, s 3 } = and [W3, s 3 ] = 0. Finally, if chiral symmetry is broken 
by terms which couple same sublattices, we must set W3 = 0. All the above conditions 
imply that the unitary matrices T belong to a group G if u = which is lowered to a group 
H at uj 7^ 0. In Table | we list the coset spaces G/H for the different classes (i) time reversal 
invariance holds with chiral symmetry 0] or without J^fl; (H) time reversal symmetry is 
broken with chiral symmetry or without [|T| ; (Hi) a magnetic field is applied in the presence 
of chiral symmetry or without it ; and finally (iv) in the presence of magnetic impurities 
with chiral symmetry or in its absence . Indeed, already at this stage we could anticipate 
the appropriate /^-functions assuming that the diffusive modes are described by a non linear 



cr-model with the proper symmetry. However, for sake of clarity, we give here a sketchy 
derivation of the a model for the simple case of real hopping [|8] . 

Within the replica trick method, we can average e~ s over disorder. The effective action 
S = So + Simp is the sum of a regular part So given by Eq.(|2|) with Uj = t (and 0y = 0) plus 
the impurity contribution S imp = —2uH 2 (^i^j) ■ Here t is the average value of the 
nearest neighbor hopping, while ut is the variance. We notice that = ^j^i, so that, 

by introducing Xf = ^ff^ff , where a and (3 is a multilabel for the different components, 
we can write Si mp = 2u 2 t 2 X^Xj a = 2u 2 t 2 J2{ij) Tr {X; L Xj). In Fourier components 

S imp = Y / E W q Tr(X q X_ q ), (3) 

V q&BZ 

where BZ means the Brillouin zone and W q = 2uH 2 (cos q x a + cos q y a) , a being the lattice 
spacing. We can decouple (J3|) by an Hubbard-Stratonovich transformation, introducing an 
auxiliary field. However, since W q = —W q+ ^^) and W q > if q is restricted to the magnetic 
Brillouin zone (MBZ), we need to introduce two auxiliary fields defined within the MBZ, 
Qo q = Ql-q and Q3q = Ql-q @ , through which 

Simp =77 E Tw~Tr [QoqQo~q + QsqQs-q] 

V q£MBZ 4KK 9 

-h E Tr [Q^-q + iQ3 q X- q -(^)} ■ ( 4 ) 

V qeMBZ 

The above expression shows that Q corresponds to smooth fluctuations of the auxiliary field, 
while Q3 to staggered fluctuations. Namely, in the long- wavelength limit, the auxiliary field 
in real space is Qr = QoR+i(~^) R Q3R, and it is not hermitean. The a model is derived by ex- 
panding the effective action for the auxiliary field Q, which is obtained after integration over 
the Grassmann variables, around the uniform saddle point in the presence of an infinitesimal 
symmetry breaking frequency. Quite generally, the saddle point Q sp may have both a t Sz 
and a r 2 Si components, £ and F respectively. However, for a ci-wave order parameter, F — 0, 



and S satisfies 1 = — J2k e'^+t, 2 ' wriere Ek = Je 2 + A|, = — 2t(cos k x a + cos k y a) is the 

k " 

band energy, and = A [cos(k x a) — cos(k y a)} is the <i-wave gap function. The above self- 
consistency equation leads to a finite density of state iVo = £/(7ra 2 t 2 ) at the chemical poten- 



tial. We can write the Q-matrix through Q R = CT^C 1 (Q sp + P R ) T R = Q(R)+CT t R C t P R T R , 
where T R belong to the coset G/H, and describe the transverse massless modes, while P R 
parametrize the longitudinal massive modes. The long wavelength action for the transverse 
modes is then obtained by integrating out the massive fluctuations and expanding in VQ(R). 
As expected, we obtain a non linear cr-model type of action 

S[Q] = J dRTr (VQ(i?)t ■ VQ(R)) 

~8^jW n / dR \ Tr [Q( R )^Q(R)} f , (5) 

where a = ^ J2k Vefc, y efc 2 +v ^ VAfc ; can be interpreted as the quasiparticle conductance in 

\ re / 

the Drude approximation, which corresponds to the spin conductance of the system II 
is a parameter related to the staggered density of states fluctuations . The last term only 
exists if chiral symmetry is present. The /3-functions for the coupling constant g = l/(27r 2 cr) 
and for the DOS in the different universality classes are known ||10|| , and are shown in Table 
|. The RG equation for c = l/(2ir 2 U) when chiral symmetry holds can be obtained by 
noticing that a + nil is not renormalized [|2|||. Therefore (3(c) = —c 2 /3(g)/(ng 2 ). In Table | 
we define T = g / (c + ng) . The scaling behavior of the DOS are extracted from the one loop 
correction to Q(R). 

According to Table |, in the zero replica limit we obtain that, if chiral symmetry is 
absent and for non magnetic impurities, the conductance vanishes, and the DOS, which is 
finite within the simplest Born approximation, is suppressed. As shown by Ref. 0, in the 
localized phase the DOS vanishes as \E\ or E 2 depending whether time reversal symmetry 
holds or not. Quite surprisingly, magnetic impurities give a derealization correction to the 
conductance, as well as a DOS enhancement. Finally, the symmetry class in the presence of 
a weak magnetic field is the same as in the standard Quantum Hall effect |1]U . 

On the contrary, if chiral symmetry is present, the conductance stays finite, or even 
increases in the presence of a magnetic field. Without magnetic field and in the absence 



of spin flip scattering, the DOS diverges approximately like p{E) ~ exp A\/— In E /E, 
with A a model dependent constant MM- By studying an on-site disordered model in the 



unitary scattering limit, a p(E) ~ 1/ (E\v? E) was found in Ref. PH . Since in the unitary 
limit, the on-site disorder effectively reduces to a random hopping, the two results should 
coincide. Indeed the leading divergence 1/E is the same, while the subleading behavior is 
not, reflecting the different way of summing up the singularities induced by disorder. 

The quasiparticle charge modes, as well as the spin modes when magnetic impurities or 
a magnetic field are present, are not described by the non linear a- model (|5]), which only 
represents the truly massless diffusion modes. Nevertheless, charge and spin conductivities, 
a c and a s , respectively, can be still evaluated through the stiffness of the corresponding 
modes, although they acquire a mass term. We have checked that this procedure gives the 
same results as the explicit evaluation of the one loop corrections to the Q-field operators 
which describe the charge and spin conductivities by introducing sources which directly 
couple to the currents ||12|| . In this way we find that the one loop corrections 5a c /a c and 
Sa s /a s coincide with 5a /a in the absence of sublattice symmetry. On the contrary, when 
sublattice symmetry holds, quasiparticle charge conductivity is reduced at one loop order 
and in the zero replica limit according to 5a c = —a c 2glns (where s is the momentum 
rescaling factor), unless spin flip scattering is present, in which case 5a c = 0. Moreover the 
spin conductance renormalizes as a in the absence of magnetic field and magnetic impurities, 
while 5a s = otherwise. Notice that while the bare a s coincides with the bare a, the bare 
a c is given by a similar expression as the bare a but without the current vertex VA/c 

We can introduce the residual quasiparticle interaction, which we assume is repulsive, in 
the spirit of the Landau Fermi liquid theory as originally done by Finkelstein in the standard 



Anderson localization Interaction generically spoils sublattice symmetry already at the 
level of the Fock's diagrams, unless for particular cases like on-site or nearest neighbor site 
interactions. Here we consider a generic case without sublattice symmetry, where just the 
three quasiparticle scattering amplitudes considered in Ref. |L4j] are important: the particle- 
hole (p-h) scattering amplitudes in the singlet, U s , and triplet, U t , channels, and the particle- 
particle (p-p) channel scattering amplitude, U c . Being the charge not a conserved quantity, 
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the singlet p-h channel has no component on the transverse massless modes, contrary to 
the p-h triplet channel. The p-p channel couples Cooper pairs of s-wave symmetry. Since 
the real part of an s-wave order parameter is in competition with the d-wave symmetry, 
only the p-p t\ channel contributes, which corresponds to fluctuations of an is pairing order 
parameter. Moreover, their contribution to a and to the DOS is exactly the same as in the 
absence of superconductivity, apart from a factor 1/2 in U c being absent the p-p ^-channel. 
Following Ref. fl4[], we find the first order interaction corrections 



2 9 nJ 3 -U, + (6) 



<r„ N "V2 2 

If time reversal invariance is broken, the U c contribution drops out. A Zeeman term re- 
moves 2/3 of the triplet contribution, while in the presence of magnetic impurities the whole 
correction vanishes. Therefore a repulsive interaction has a delocalizing effect which com- 
petes with quantum interference corrections. Moreover, both the p-h triplet channel and 



the p-p zs-channel [HJ enhance the DOS. However, being the interaction correction pro- 
portional to the small DOS induced by disorder, it is likely negligible, although it might 
become dominant at very low temperatures if the enhancement of the p-h triplet scattering 



amplitude takes place WM, while the p-p one is suppressed at higher energy by the Cooper 



phenomenon []14| , [T5| . For large values of the triplet amplitude Eq. (|6]) is no more appropriate 
and a RPA resummation of the dynamical amplitude should be carried out leading to the 
same expressions for <5cr/cr and 6N/N Q as in the absence of superconductivity fll4f . 

Finally, let us compare the above theoretical predictions with the experimental behavior 
of the superconducting cuprates. In the realistic situation in which time reversal symmetry 
is unbroken and chiral symmetry is absent, the quantum interference corrections should 
appear when In r^/r ~ (vp + v\)/vpVA, where r = 1/2S is the elastic relaxation time 
induced by disorder, the inelastic one, Vp and v& the velocities parallel and perpendicular 
to the nodal directions. Due to the still controversial situation for the inelastic relaxation 
time in BSCO, let us estimate the crossover temperature T* below which localization effects 



should appear in YBCO. On the basis of the microwave conductivity data of Ref. |L6j for 



optimally doped YBCO, one finds T* ~ 1 K if oc T 42 is assumed [JlB] , or ~ 0.3 K 



for oc T~ 3 |[L7| . Moreover, T* might be further reduced either by interaction and nesting 
effects, as shown previously, or by a finite correlation range of the impurity potential [|T8|,|T|| . 
Experimentally, optimally doped YBCO and BSCO thermal conductivities [0,0 do not 
show any localization suppression down to 0.1 K, which is however not so small as compared 
to our estimates for T* to exclude that localization finally takes place. Indeed in underdoped 
YBCO there is evidence of a vanishing residual quasiparticle conductivity, as extracted by 
heat transport . 



In summary, we have derived the general form of the non linear a-model, which describes 
the quasiparticle diffusive modes in a c?-wave supercondutor, in the presence of a random 
hopping. Away from half-filling, the random hopping plays a similar role as a random on-site 
potential. However, at half-filling, the model posseses an additional chiral symmetry and, 
as a result, it describes a critical theory with a finite spin conductance. We believe that the 
results are of more general validity than the peculiar choice of the disorder would suggest. 
Indeed the model is also applicable in the case of a strong disordered impurity potential close 
to the unitary scattering limit, which effectively reduces to a random hopping. Incidentally, 
this is believed to be just the case for the High T c superconductors. We have also calculated 
the localization corrections to the quasiparticle charge conductivity, namely of the optical 
conductivity at small frequency, and find that it should vanish at sufficiently low temperature 
in all cases but in the presence of magnetic impurities. 

Finally, we have included in our analysis a residual quasiparticle interaction, which we 
have shown has a delocalizing effect mainly due to the quasiparticle spin diffusive modes. The 
results suggest that these systems, although being singlet superconductors, might support 
consistent spin fluctuations. 

Note added While writing this paper, we became aware that Jeng, Ludwig, Senthil and 
Chamon |22[ have studied the interaction effects in the same model within the Keldisch 



technique, obtaining results similar to ours. 
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TABLES 

TABLE I. Coset spaces and (3 functions for the coupling g and for the DOS N in the different 
universality classes. T is the time reversal invariance. 





Coset space 


1 y 


/3jv 


Yes Chiral, Yes f 


U(4n)xU(4n)/U(4n) 


8na 2 


(T/4 - 8n)o 


Yes Chiral No f 

_i_ v^kj v_- 1111 capi - 11 v_y -J- 


JJ(4ri)/0(4ri) 


Ana* 


f-i + r/4 - 471)0 


Yes Chiral, magnetic field 


0(4n)/0(2n)xO(2n) 


(2n-l)g 2 


-2ng 


Yes Chiral, spin flip 


U(2n)/U(n)xU(n) 


ng 2 


-ng 


No Chiral, Yes f 


Sp(2n)xSp(2n)/Sp(2n) 


2{2n+l)g 2 


(-l-4n)g 


No Chiral, No f 


Sp(2n)/U(2n) 


(2n + l)g 2 


(-l-2n)g 


No Chiral, magnetic field 


U(2n)/U(n)xU(n) 


ng 2 


-ng 


No Chiral, spin flip 


0(2n)/U(n) 


(n - 1) 5 2 


(1 - n)# 
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